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Spin chain in magnetic field: limitations of the large -N mean-field theory 
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Motivated by the recent success in describing the spin and orbital spectrum of a spin-orbital 
chain using a large-A mean-field approximation [J, we apply the same formalism to the case of 
a spin chain in the external magnetic field. It occurs that in this case, which corresponds to 
N = 2 in the approximation, the large-A mean-field theory cannot qualitatively reproduce the spin 
excitation spectra at high magnetic fields, which polarize more than 50% of the spins in the magnetic 
ground state. This, rather counterintuitively, shows that the physics of a spin chain can under some 
circumstances be regarded as more complex than the physics of a spin-orbital chain. 
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Introduction Recently a number of studies discussed 
the collective excitations in a spin-orbital chain [mans 
H 0. Most of them concentrated around a novel phe¬ 
nomenon called spin-orbital separation which is present 
when a very strong external crystal field fully polarizes 
the orbital sector of the ground state mm- Although 
this phenomenon seemed to be completely at odds with 
the physics present in an SU(4)-symmetric spin-orbital 
chain (i.e. without external crystal field) [6], a very re¬ 
cent paper discusses how to unify these two seemingly 
different limits [1 . It occurs that a large-A mean-field 
theory mm surprisingly well describes the spin and or¬ 
bital spectra for any value of the crystal field and thus 
explains the striking evolution of the spin and orbital 
spectra with increasing external crystal field [1 . 

As a result of this recent success of the large-A mean- 
field theory the following question arises: could such a 
theory be equally successful in explaining the behavior of 
collective spin excitations in a spin chain that is subject 
to external magnetic field? While this might look like as 
an old problem, which should have been solved long time 
ago, to the best of our knowledge, there exists no precise 
answer to this question in the literature. 

Definition of the problem Let us now be more spe¬ 
cific. First, we define the following Hamiltonian which 
describes the problem of a spin chain subject to the ex¬ 
ternal magnetic field H z 

H=j'52(SiS j + \)+H z Y, s *- (!) 

(U) 1 

Here J is the energy scale of the superexchange inter¬ 
actions between SU(2)-invariant spin 5 = 1/2 operators 
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(S), (ij) represents a nearest-neighbor spin pair, and H z 
is the magnetic field strength. 

Second, we define the transverse dynamical spin struc¬ 
ture factor which is a good proxy for probing the nature 
of the collective spin excitations: 

s (»■">- 1 (2 > 

Here \^) is the ground state of H with energy S% = 

JT e iq i S? /\[L is the Fourier transform of the local spin 
operator, and L is the number of lattice sites. 

In what follows we calculate the dynamical spin struc¬ 
ture factor in Eq. Q using two distinct methods: (i) 
the numerically exact combined cluster perturbation the¬ 
ory (CPT) and exact diagonalization (ED) method, and 
(ii) the approximate analytical large-A mean-field the¬ 
ory. We compare our analytical results with the numeri¬ 
cal calculations, which have already been extensively dis¬ 
cussed in the literature 

Numerical results The numerical method employed in 
the current study, CPT+ED, is a quantum cluster ap¬ 
proach m urn which complements the finite-size ED 
simulations and therefore allows for a better visualiza¬ 
tion of the fine spectral details. The spin dynamical 
structure factor, Eq. (J2|, calculated with this method is 
shown in Fig. [l] for three different values of the external 
magnetic field H z . When H z = 0, the ground state has 
short range antiferromagnetic order and the spectrum is 
well-known ia uni nu ns]: it is mostly spanned by a 
two-spinon continuum and has zero modes at q = 0 and 
q = 7r, cf. Fig. [lja). For a finite value of H z the ferro¬ 
magnetic domains start to appear in the ground state and 
the positions of the zero modes shift. For instance when 
H z ~ 1.58J half of the spins in the chain are polarized 
and the spectrum has zero modes at q = n/2 and q = 7r, 
cf. Fig. 0b). Nevertheless, the spectrum is still spanned 
by a continuum of fractional excitations. This, however, 
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FIG. 1: Numerical results Spin dynamical structure factors computed by CPT+ED for a spin chain under a magnetic field 
H z : (a) H z — 0 with no spin polarization (JT Sf /L m 0); (b) H z ~ 0.79 H c fi with half polarized spins (JT Sf/L m 1/4); 
(c) H z — H z r — 2.0 J with fully polarized spins (JT Sf /L = 1/2). The ED spectra (broadened with a 0.25J Lorentzian) are 
computed on an L = 24 site lattice. For H z < H z r [panels (a) and (b)], the spin spectra show fractionalized excitations with 
broad energy continua. When H z > i/J r , the (ferromagnetic) ground state is fully spin-polarized; the excitations are no longer 
fractional, and the spectrum exhibits only a sharp, single-magnon mode [panel (c)]. 


stays in contrast with the fully polarized (i.e. ferromag¬ 
netic) ground state which occurs at H z = 2 J = H^ r . The 
spin dynamical structure factor is then no longer spanned 
by a two-spinon continuum, cf. Fig. [IJc), and instead a 
single (magnon) branch arises. Further increasing H z 
above H% r leads to a gapped magnon branch. 

Large-N mean-field theory Following similar steps as 
described in detail in Refs. pans] we first map the spin 
model described by Eq. 0 onto a fermionic model m- 
By performing a large- N mean-field decoupling and solv¬ 
ing the self-consistent mean-field equations, we obtain 
the following Hamiltonian 

^MF = + Skifljki ), (3) 

k 

where the fermionic bands are / 1 = 

—2 J cos(Sk) cos k/'K =p Jsin(25fc)/7r (cf. Fig. [2| with 
5k = arcsin[i7 2 7r/(2J)]/2, subject to the constraint 

We note, first, that in this approach H z < 2 J/n, which 
means 5k < tt/4 and a maximum Fermi momentum 5k + 
kp = 37r/4 under an applied magnetic field (since kp = 
7r/2 for the fermionic mean-field theory at H z = 0 [T6]). 
Second, above H z > 2J /7r the self-consistent mean-field 
solution effectively breaks down. 

Next, since the maximum available Fermi momentum, 
37r/4, corresponds to just half of the spins being polarized 
in the ground state, this is the maximum value of polar¬ 
ization available in this approach. As this means that we 
can never fully fill one of the spin bands, i.e. reach the 
5k ~F kp =7r momentum, the fermionic mean-field the¬ 
ory is not able to describe the fully spin-polarized ground 
state. 

In this fermionic mean-field picture, the compact sup¬ 
port of the spin spectrum can be calculated by S(q , uj) = 

^keFS,q+k^FS ~ £ q+k,t + £ ki) + J2keFS,q+k(£FS 
£ q +k,i + c f- Fig- ^1 The evolution of the spin spec¬ 
trum as a function of H z is shown in Fig. [3|a)-(b). The 
shift of the zero-energy modes with the magnetic field 




FIG. 2: Large-N mean-field theory Evolution of the mean- 
field fermionic bands as a function of the magnetic field at 
H z — 0 (top panel) and H z m 2J/7T (bottom panel). The 
collective spin excitations in the mean-field picture correspond 
to ‘particle-hole’ excitations of the fermions across the Fermi 
surface (denoted by the dotted horizontal lines). The energies 
of the up-spin and down-spin fermionic bands are separated 
by H z , and the allowed ‘particle-hole’ excitations change with 
the magnetic held accordingly. The thick arrows point to the 
allowed zero-energy spin excitations. 
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FIG. 3: Analytical results Compact support of the spin spectra for a spin chain in a magnetic field H z computed by the large-iV 
mean-held theory [(a)-(b)]: (a) When H z = 0, the ground state exhibits (short-range) antiferromagnetic correlations without 
any spin polarization, (b) When H z = 2J/7T in the mean-held calculation, i.e. when 50% of the spins are polarized in the 
ground state; the lighter (darker) part in the spin spectrum refers to spin-hip excitations created by spin raising (lowering) 
operators, (c) The exact spin wave dispersion oo q — J(l + cosg) calculated by linear spin wave theory at H z — H^ r \ the large -N 
mean-held theory is not valid in this regime. 


in the spin spectrum follows from the splitting of the 
fermionic bands in the magnetic held, cf. Fig. [2] 

Conclusions Let us hrst emphasize that the fermionic 
large-7V mean-held approximation is valid only when 
H z < 2J/7T (i.e. when the ground state is less than half- 
polarized). In this case there exists a qualitative agree¬ 
ment between the analytically calculated compact sup¬ 
port and the numerically obtained spin spectrum. Nev¬ 
ertheless, the quantitative differences between the mean- 
held and the CPT+ED results are substantial (the band¬ 
width of the spin excitations is ca. 2.5 smaller in the an¬ 
alytical approach), making the fermionic large -N mean- 
held description less appropriate compared to that in the 
spin-orbital model pQ. 

A different situation exists when H z > 2 J/tt (i.e. when 
the ground state is more than half-polarized). Here the 
mean-held approximation breaks down. This means that 
long before the critical held is reached (which fully po¬ 
larizes the ground state) there exists no simple analytical 
approach that can be used to calculate the spin spectrum 
(for more complex approaches, cf. Ref. [9} HZ]). It is then 
only when H z > H^ r that another simple analytical ap¬ 
proach becomes valid - the well-known linear spin-wave 
approximation which works well for spin ground states 
with long range magnetic order mm- The analyti¬ 
cally evaluated spin spectrum, S(q 1 oo) = S(oo — oo q ) with 
ou q = J( 1 + cosg) being the magnon dispersion when 
H z = H% r [cf. Fig. [3^c)] , is then equal to the numerically 
calculated spin dynamical structure factor [cf. Fig. [ljc)]. 
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